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摘 Ze. LL 波 段 太 阳 射 电 爆 发 是 导航 系统 不 稳定 性 的 潜在 影响 因素 ， 通 过 工 波 段 内 精密 太阳 射电 流量 的 
监测 可 以 实时 发 现 太 阳 射 电 爆 发 干扰 导航 事件 ， 为 此 云南 天 文 台 本 部 拟 建立 一 个 L 波段 多 频 点 太阳 射电 监 
测 系 统 。 而 对 无 线 电 环 境 的 有 效 评估 对 于 该 监测 系统 观测 数据 的 稳定 获取 至 关 重 要 。 本 文 介绍 了 该 监测 平 
台 的 无 线 电 监测 准备 研究 ,通过 对 云南 天 文 台 本 部 凤凰 山地 区 的 工 波段 无 线 电 环境 进行 了 100 小 时 的 测试 。 
本 文 提出 一 种 基于 Simple Thresholding 算法 与 CUSUM 算法 的 改进 辣 值 的 算法 ， 送 选 出 介 于 北斗 Bl. B2. 
B3 频 点 ，GPS Ll. L2 频 点 之 间 7 个 5MHz 无 线 电 干扰 较 少 的 无 线 电 通 带 ， 分 别 为 1551MHz~1555 MHz, 
1596MHz~1600MHz、 1161MHz~1165MHz、 1221 MHz-1225MHz., 1246MHz-1250MHz, 1291MHz-1295MHz., 
1231MHz-1235MHz, i$: 3E 2) 30] 7g 98.32996., 98.301906. 98.315906, 98.33596., 98.224% 97.650% 98.260906, 
均 符 合 太阳 观测 需求 。 为 下 一 步 接收 机 的 设计 和 信号 处 理 提 供 了 依据 。 
关键 词 : LRR; 监测 系统 ; 无 线 电 环境 ; BMB. 洁净 率 
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太阳 活动 是 空间 天 和 气 的 源头 ， 日 铣 物 质 抛 射 (coronal mass ejection, CME) pao 
REDE. SB HEOUUU Sc RR T AS F8] ELSE REATARD, uda HE. Ries, Nt—c:1b. BESSER TA 
从 太阳 上 抛射 出 来 ， 其 密度 随 着 远离 太阳 的 距离 而 降低 随 着 CME 物质 远离 太阳 ， Rd 电 频 率 
也 在 随 之 降低 ， 为 此 太阳 物理 界 一 致 认 为 地 基 太 阳 射 电 望 远 镜 所 能 观测 到 的 无 线 电 爆 发 频率 越 低 ，CME 物 
质 离 地 球 就 越 近 。 那 么 在 不 同 频 率 上 的 太阳 射电 进行 观测 ， 必 然 可 以 反 演 这 个 日 时 高 度 上 的 动态 信息 中 1。 
同时 ， 在 长 期 的 太阳 射电 流量 监测 中 发 现 其 不 同 频段 的 流量 变化 与 太阳 活动 周 高 度 相 关 ， 其 是 重要 的 
空间 天 气 观测 指标 。 但 另 一 方面 ， 射 电 的 观测 频段 与 民用 无 线 电 使 用 频段 重合 ， 这 样 极 易 容 易 造 成 对 观测 
数据 的 干扰 ， 造 成 误 判 。 在 射电 望远镜 建设 之 前 的 无 线 电 环境 测试 与 评估 尤为 重要 ， 与 脉冲 星 、 活 跃 星系 
核 (active galactic nucleus，AGN) 光 变 等 需要 宽带 观测 不 同上 1。 由 于 高 时 间 分 辨 率 要 求 ， 太 阳 射 电 精密 流量 
观测 时 对 一 些 侦 发 性 的 干扰 并 不 敏感 ,这 个 短 时 “ 坏 点 ”可 以 通过 后 期 与 太阳 光学 等 数据 进行 比较 排除 爆发 可 


eu 
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L 波段 太阳 射电 爆发 是 导航 通信 的 潜在 干扰 因素 , 与 其 他 空间 天 气 事件 不 同 的 是 , 太阳 爆发 的 强 信和 号 以 
光速 传播 ， 是 作用 最 快 的 空间 天 气 事 件 ， 且 从 理论 上 来 说 整个 迎 日 面 的 导航 接收 机 均 有 可 能 受到 干扰 9， 
董 亮 、 黄 文 耿 、 闫 小 娟 等 人 提出 了 采用 L 波段 多 通道 精密 太阳 射电 望远镜 同步 监测 的 办 法 ， 实 时 监测 工 波 
段 太 阳 射 电流 量 ， 在 爆发 早期 发 出 预警 信号 5 。 

但 是 工 波 段 存在 大 量 无 线 电 干扰 , 在 建设 望远镜 前 本 文 需要 评估 MHz 级 的 观测 带 内 无 线 电 干扰 信号 的 
分 布 规律 ， 判 断 可 观测 太阳 时 间 ， 对 于 长 时 间 被 干扰 占用 的 信道 进行 排查 。 以 防 发 生 误 报 等 情况 ， 针 对 即 
将 在 云南 天 文 台 凤 凰 山本 部 建设 的 工 波 段 的 精密 流量 太阳 射电 望远镜 (针对 中 高 日 冕 进行 研究 和 太阳 射电 爆 
发 干扰 导航 通信 预警 ， 并 作为 廊坊 太阳 射电 望远镜 的 备份 )， 本 文 介 绍 了 对 台 址 的 无 线 电 测试 结果 和 观测 频 
段 遂 选 方法 以 及 遂 选 出 介 于 北斗 B1、B2、B3 频 点 ，GPS L1. L2 频 点 之 间 7 个 5MHz 无 线 电 干 扰 较 少 的 
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无 线 电 通 带 。 


1 系统 组 成 
1.1 系统 组 成 

测试 系统 由 HL050 天 线 、 低 噪声 放大 器 、 稳 压 电 源 、FSU256 频谱 仪 和 计算 机 组 成 ， 测 试 系统 框图 
如 图 1 所 示 。 
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图 1 测试 系统 
Fig.1 Test system 

系统 框图 中 ,由 HL050 天 线 将 采集 的 信号 通过 低 噪 声 放 大 器 放大 信号 ， 再 传输 到 FSU26 频谱 仪 转化 成 
频谱 ， 并 把 数据 传输 给 计算 机 进行 显示 和 存储 。 其 中 稳 压 电源 为 放大 器 提供 稳定 电压 , 保证 数据 的 准确 性 ， 
计算 机 通过 LabView 图 形 化 编程 实现 控制 频段 范围 和 极 化 方式 ， 主 要 分 为 两 个 部 分 : HL050 RRE E 
和 频谱 仪 控制 程序 "”。 
1.2 测试 设备 及 参数 

天 线 采 用 的 是 HL050 对 数 周期 天 线 ， 参 数 见 表 1，HL050 天 线 是 一 种 应 用 广泛 的 天 线 ， 其 优点 是 结构 
简单 、 频 带宽 、 功 率 容 量 大 、 调 整 与 使 用 方便 "。 

频谱 仪 采用 德国 罗 德 施 瓦 茨 公司 的 FSU26， 工 作 频 率 20Hz~26.5GHz， 平 均 噪声 电 平 为 -138dBm2， 参 
数 见 表 1。 

计算 机 利用 LabView 软件 图 形 化 编程 实现 ， 在 LabView 软件 图 形 化 编程 中 ， 可 以 调节 测试 的 频段 起 始 
频率 、 极 化 方式 、 分 辩 率 带宽 、 积 分 时 间 及 执行 次 数 ， 并 将 采集 到 的 信号 保存 在 预 设 地 址 中 ， 便 于 对 后 续 
数据 的 处 理 分 析 。 
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Tab.1 Test equipment and parameters”! 


Equipment Model Frequency range Technical specifications 
R&S logarithmic periodic antenna HL050 0.85MHz-26GHz Standing Wavebie: <2.5 gain /dBi: 8-8.5 
B&Z low noise amplifier BZ-0218-201030-202525 2GHz-18GHz Standing Wavebie: <2.5 gain /dBi: 30 
R&S FSU26 spectrometer RS-HL050 20Hz-26.5GHz Phase noise: 137dBc /Hz 


1.3. 监测 过 程 
测试 地 点 为 云南 省 昆明 市 中 国 科 学 院 云南 天 文 台 , 主要 对 天 文 台 凤凰 山本 部 天 文 台 台 址 ( 25°1'48.538"N , 
102°48'14.180"E ) 东 、 西 、 南 、 北 4 个 方向 工 波段 的 测试 ， 系 统 频率 范围 在 1GHz~2GHz， 频 谱 仪 扫 频 宽度 
设置 为 1000MHz。 频 谱 仪 显示 带宽 和 分 辨 率 带 宽 设 置 为 300KHz， 扫 描 时 间 30s。4 个 方向 都 测试 了 4 次 且 
时 间 超 过 8 小 时 ， 共 16 次 测试 ， 每 组 数据 均 超 过 800 组 。 有 具体 测试 日 期 及 其 组 数 如 表 2 所 示 。 
表 2 测试 日 期 及 组 数 (2022 Æ) 
Tab.2 Date of test and number of groups (in 2022) 


Test the orientation First Second Third fourth 
East 04/04 831 group 04/13 1306 groups 04/24 1175 groups 04/25 1088 groups 
South 04/03 820 groups 04/14 1368 groups 04/23 924 groups 04/26 1117 group 


West 04/07 1062 groups 04/15 882 groups 04/21 988 groups 04/27 2219 Group 
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North 04/11 1074 groups 04/12 1270 groups 04/22 1137 groups 04/28 804 groups 


2 数据 处 理 


本 节 2 NOB 28 s BE REE Cg E RARR AE E E E LTEN A E o 
2.1 数据 介 

文件 的 格式 ， 单 位 为 qBi ， 缩 小 数据 的 绝对 数值 ， 方 便 记 录 。 由 于 整 
个 测试 过 程 设 置 频谱 仪 输出 两 条 线 ， 所 以 每 一 组 检测 数据 里 面 分 为 两 列 ， 如 图 2 如 示 。 第 一 列 为 所 采集 交 
据 的 最 大 保持 值 ， 第 二 列 为 所 扫描 时 间 内 的 平均 值 ， 为 每 次 检测 时 判断 测试 系统 的 稳定 效果 。 且 数据 储存 
为 非 线性 的 数据 格式 ， 这 是 一 种 对 数 的 关系 ,底数 为 10， 数 学 关系 为 Jog10(P /P), IF P/P 表示 功率 比 。 
故 后 续 数据 处 理 需要 对 数据 进行 数据 的 线性 转化 0627。 
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frequency (MHz) 
图 2 测试 系统 实测 图 


Fig.2 Actual measurement diagram of the test system. 


2.2 WERENT 9e 
针对 即将 在 云南 天 文 台 凤凰 山本 部 建设 的 工 RRR AREA E, Ade 
B1(1559.052 MHz~1591.788 MHz), B2(1166.220MHz-1217.370 MHz). B3(1250.618 MHz~1286.432 MHz) 频 


干扰 较 少 的 无 线 电 通 带 ， my GPS 中 L1 频 点 在 北斗 Bl 频 点 的 范围 中 ， 且 B2 与 12 5 点 之 间 相差 没 
有 超过 10MHz, 无 法 遂 选 两 处 相对 洁净 的 5MHz 的 通 带 ,B2 右边 送 选 的 通 带 与 L2 左边 遂 选 的 通 带 相 重 合 ， 
一 共 7 o 道 〈 经 过 无 线 电 环境 监测 带 内 无 干扰 信号 )， 使 得 GPS、 北斗 的 所 有 通信 频 点 都 穿插 在 监测 频 点 
之 间 。 

本 次 测试 的 目的 是 从 工 波段 中 六 选 出 洁净 率 〔 无 线 电 干扰 的 时 间 占 总 监测 时 间 的 百分比 高 于 95% 的 
可 用 信道 。 遂 选 方式 则 是 通过 绘制 出 的 频谱 瀑布 图 中 能 直观 的 选 出 较为 干净 的 信道 ， 通 过 遍历 观测 数据 超 
过 阔 值 的 占 比 来 确定 频谱 洁净 率 , 当 遂 选 出 的 信道 的 洁净 率 都 大 于 95% 则 说 明 该 信道 可 用 。 处 理 流程 如 下 : 
入 一 步 ， 先 选取 一 个 短 时 “宁静 ”的 观测 频段 作为 闵 值 计算 数据 样本 ; 

第 二 步 ， 在 上 述 数 据 样本 中 通过 本 文 提出 的 闵 值 算法 计算 ， 得 到 病 值 ; 

币 三 步 ， 以 此 闵 值 再 裔 历 其 余 观 测 时 间 段 数据 ， 对 比 赚 值 ， 超 过 羡 值 时 间 和 低 于 阔 值 时 间 进 行 对 比 ， 
求 取得 到 该 时 间 段 的 干扰 信号 所 占 比例 ; 

第 四 步 ， 利 用 算法 计算 得 到 的 阔 值 对 本 文 遂 选 的 7 个 5MHz 无 线 电 干扰 较 少 的 无 线 电 通 带 进行 评估 。 
22. BB 5 

闵 值 分 析 法 因 其 实现 简单 、 检 测 结果 精度 较 高 而 被 广泛 应 用 。 代表 性 阔 值 分 析 法 有 CUSUM (cumulative 
sum) 法 中 和 Simple Thresholding 法 请， 在 检测 某 行 ( 列 ) 的 观测 数据 时 ，CUSUM 算法 是 通过 估算 累积 样本 的 
方差 与 平均 值得 到 阔 值 , Simple Thresholding 算法 是 使 用 该 行 ( 列 ) 的 中 位 数 作为 其 阔 值 的。 

首先 本 文 监 测 过 来 的 数据 取 对 数 之 后 的 数据 ， 取 对 数 之 后 不 会 改变 数据 的 性 质 和 相关 关系 ， 且 压缩 了 
变量 的 尺度 ， 数 据 更 加 平稳 ， 也 消 弱 了 模型 的 共 线 性 、 异 方差 性 等 ， 所 以 计算 闷 值 时 需要 将 数据 线性 化 ， 
如 式 (1) 所 示 。 
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(1) 


样本 的 方差 是 表示 RFI 的 统计 异常 值 的 重要 估计 值 。 由 于 绝对 中 位 差 (median absolute deviation; MAD) 


和 中 位 数 是 消除 异常 值 的 最 有 效 估计 ， 因 此 本 文 提 出 了 一 种 基 


于 Simple Thresholding 算法 与 CUSUM 算法 


的 改进 阔 值 算法 对 台 址 的 L 波段 无 线 电 测试 结果 和 观测 频段 i 


行 闭 选 ， 如 果 所 选取 信道 闪 值 水 平 超过 7h ， 


WL RFI. Th 如 式 (2) 所 示 。 
Th-X — a4MAD 


(2) 


式 (3) 中 , Th 为 CUSUM SEHE BISRZRTEBIB, X 为 线性 数据 的 均值 ,其 中 0.0 «ac < 5.0, 一 般 默 认 ac =3.0; 


MAD - 1.4286» med ,,,... ( 


H. med — 0.5(x, + x, 


m+ 


x, - med(x,)]) 


(3) 


), n-2m; med-x,,,» n-2m«l. Ñ}, MAD 是 Simple Thresholding 算法 


中 使 用 的 方差 估计 值 ， 由 等 式 (3) 定 义 ， 其 中 中 位 数 是 原始 总 体 数 组 x 的 中 位 数 。 然 后 用 此 中 位 数 减 去 原始 
样本 每 个 元 素 ， 创 建 与 原始 元 素 大 小 相同 的 新 修改 样本 。 然 后 计算 这 个 新 样本 的 中 位 数 ， 并 用 恒定 的 比例 


因子 1.4286 相 乘 号 ， 以 使 该 估计 与 预期 的 高 斯 分 布 一 致 。 


得 到 线性 计算 的 闵 值 Th 之 后 ， 需 要 将 阀 值 取 对 数 转换 在 监 


监测 的 数据 中 去 计算 通 带 的 洁净 率 ， 如 式 (4) 


所 示 ， 得 到 最 后 的 阀 值 结果 7h* 。 整 个 阀 值 算 法 框图 以 及 算法 伪 代 码 如 图 


Begin ) 
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z Read in the data 7 
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Y 
MAD - 1.4286 x med,.....( 


X% = med(x,)|) 


Th* 210? *log10(Th) 


3 所 示 。 


(4) 


Algorithm 1: Threshold algorithm 
Data: current data T;, Linear data z;,Linear threshold T'h,Final 
threshold 
Result: Channel threshold Th* 
1 R= (0,5...) 


Es 2 rj = (i, 25..4]); 
X, (1*10 3)10 10 AT s P 2.2: 
4 Ari t []; 


5 while zi = (1 1073)1087:/10 and AT; 4 0 do 

6 z; = med(zi)i=(14n)/2 TQ) € {21 < T2 €. < En}; 

7 Arz; = (|ri — z+|} Az; € (Azi, Azo, ..., Arn}; 

8 | X (nini. z5)/n Find the linear mean; 

9 MAD = umed(Az;) Ji = 1.4286 MAD is an estimate of variance 
used in the Simple Thresholding algorithm; 


i Y ic 10 Th = X -aMAD Th is the linear threshold of the CUSUM 
一 algorithm; 
MEX SAMAD i | if 0.0 <a < 5.0 then 
y 12 & = 3.0; 
[m -10^*1og10CTh) 13 Th* — Th; 
14 end 
Y 15 Th* = 107? « log10(Th) The algorithm finds the final threshold; 
A Output Th* 7 16 end 


End 


(a) 


图 3 闵 值 算法 框 


名 | 


(b) 


法 伪 代 码 .(a) RERE; (b) 闵 值 算法 伪 代 码 . 


Fig.3 Threshold algorithm block diagram and algorithm pseudocode .(a) Threshold algorithm block diagram; (b) Threshold algorithm pseudocode. 


2. 2. 


2 信道 选取 与 数据 比较 
第 一 步 , 先 选 取 一 个 短 时 “ 宁 前 
最 多 的 一 次 作为 信道 的 初 选 ， 西 


向 测试 在 第 四 次 04 月 27 


”的 观测 频段 作为 阀 值 计算 数据 样本 , 初 选 是 从 表 2 中 选取 此 次 检测 组 数 


日 时 共 采 集 了 2219 组 数据 ， 故 此 次 选取 短 时 “ 宁 
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静 ” 的 观测 频段 是 在 西向 第 四 次 数据 所 画 的 频谱 瀑布 图 中 选取 ， 频 谱 图 如 图 4~ 图 6 所 示 ; 

第 二 步 ， 在 上 述 数 据 样本 中 通过 上 述 的 闵 值 算法 计算 得 到 冰 值 ; 

第 三 步 ， 以 此 闵 值 再 吉 历 其 余 观 测 时 间 段 数据 ， 对 比 闵 值 ， 超 过 阔 值 时 间 和 低 于 羡 值 时 间 进 行 对 比 ， 
到 该 时 间 段 的 干扰 信号 所 占 比例 ， 对 比 结果 如 表 3 所 示 ; 


求 取 


得 
第 四 步 ， 利 用 算法 计算 得 到 的 阔 值 对 本 文 交 选 的 7 个 SMHz 无 线 电 干扰 较 少 的 无 线 电 通 带 进行 评估 。 


L-band test data 


Signal strength 
I] 
[9] 


1800 ^ —— ~ 2000 
Spectral channels 2000 2500 Number of samples 


图 4 L 波段 频谱 瀑布 图 
Fig.4 L-band spectral waterfall 
4 中 横 坐 标 为 工 波段 1GHz~2GHz ， 纵 坐标 为 数据 采样 的 组 数 ， 竖 轴 为 采样 的 的 信号 强度 ， 因 为 能 
较为 直观 的 从 图 中 看 出 本 文 所 需 频段 的 初步 干扰 情况 ， 故 此 次 工 波 段 的 所 有 测试 均 采 用 此 种 绘制 方法 。 


L-band test data 


-20 ~ 


-30 - 


-0r- 


-50 - 


Signal strength 


-90 | | | | 
1000 1200 1400 1600 1800 2000 
Spectral channels 


图 5L 波段 频谱 右 视图 
Fig.5 L-band spectrum right view 
图 5 73 L 波段 频谱 瀑布 图 的 右 视 图 ， 图 中 能 清晰 看 到 所 有 测试 组 数 的 信号 强度 ， 且 能 直观 看 到 本 文 所 
需 的 频段 的 干扰 情况 ， 能 为 本 文 寻找 短 时 “宁静 ”的 观测 频段 确定 的 初步 第 选 工 作 。 图 6 为 L 波段 频谱 瀑布 
图 的 下 视图 ， 图 中 标记 出 北斗 B1. B2. B3 频 点 及 GPS L1. L2 频 点 位 置 ， 为 后 续 通 带 的 遂 选 起 辅助 作用 。 
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L-band test data 


Ta 
28 


Spectral channels 


D J 
四 pS 
BEGES 8 


E 


: 


500 1000 1500 2000 2500 
Number of samples 


图 6L 波段 频谱 左 视图 
Fig.6 L-band spectrum left view 
通过 阔 值 算法 对 数据 闵 值 的 确定 ， 实 现 对 整个 L 波段 监测 数据 六 选 的 多 个 5MHz 无 线 电 干扰 较 少 的 无 
线 电 通 带 中 确定 最 终 效 选 的 介 于 北斗 星 及 GPS 频 点 附近 的 7 个 无 线 电 通 带 如 图 8 所 示 ， 具 体 的 信道 选取 的 
流程 框图 如 图 7 所 示 。 本 文 初步 六 选 了 13 个 通 带 ,但 经 过 最 终 信 道 选 取 之 后 确定 了 7 个 5MHz 无 线 电 干扰 
较 少 的 无 线 电 通 带 。 


Input the data 
Choose multiple sets of quiet bands 上 ----- » The groupings are calculated sequentially 


>| Linear transformation of the data — 


! 


Give the threshold 


Whether the cleanliness rate is qualified? 


Yes 


Continue to the next set of decisions 


This band is rejected 


No 


Outputs this frequency band 


图 7 送 选 信道 流程 框图 
Fig.7 Block diagram of the selection channel flow 

本 次 测试 主要 是 为 了 遂 选 出 介 于 北斗 Bl. B2. B3 频 点 ，GPS L1. L2 频 点 之 间 7 个 5MHz 无 线 电 干 

扰 较 少 的 无 线 电 通 带 ， 首 先 从 图 4~ 图 6 中 L 波段 频谱 图 中 初 选 出 5MHz 无 线 电 干扰 较 少 的 无 线 电 通 带 

1161MHz-1165MHz 频段 作为 初 选 闪 值 确定 的 可 用 信道 ,再 根据 上 述 算 法 对 阔 值 的 确定 后 ,通过 闵 值 遍历 对 

比 程序 将 考察 信道 的 洁净 率 进 行 计 算 ， 得 到 各 组 监测 数据 考察 信道 的 洁净 率 如 表 3 所 示 。 从 表 3 中 可 以 看 

出 在 总 的 16 次 工 波 段 的 测试 中 所 遂 选 的 SMHz 无 线 电 干扰 较 少 的 无 线 电 通 带 1161MHz~1165MHz 频段 的 信 
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道 洁 净 率 都 在 95% 以 上 ， 符 合 太 阳 观 测 需 求 ， 所 以 本 文选 用 此 阔 值 为 所 送 选 的 7 个 5MHz 无 线 电 干扰 较 少 
的 无 线 电 通 带 的 阔 值 。 
表 3 1161MHz~1165MHz 频段 洁净 率 
Tab.3 Cleanliness rates in the 1161MHz-1165MHz band 


Test the orientation First Second Third fourth 
East 98.580% 98.285% 98.366% 98.106% 
South 97.829% 98.509% 98.051% 98.388% 
West 98.418% 98.413% 98.744% 98.215% 
North 98.585% 98.472% 98.205% 97.869% 


本 文 结合 图 5~ 图 7 以 及 表 3 初步 遂 选 介 于 北斗 B1. B2. B3 频 点 以 及 GPS L1, L2 频 点 之 间 7 个 5MHz 
无 线 电 干扰 较 少 的 无 线 电 通 带 如 图 8 所 示 ， 图 8 73 L 波段 频谱 瀑布 图 中 所 遂 选 的 通 带 情况 ， 从 图 中 看 到 本 
文 所 选 的 通 带 具体 的 干扰 情况 : 

图 中 B2a(1161MHz~1165MHz) 和 B2b/L2a(1221MHz~1225MHz) 为 B2 频 点 附近 通 带 ; 

B2b/L2a(1221MHz~1225MHz) 和 L2b(1231MHz~1235MHz) 为 L2 频 点 附近 通 带 ; 

B3a(1246MHz~1250MHz) 和 B3b(1291MHz~1295MHz) 为 B3 频 点 附近 通 带 ; 

Lla(1551MHz~1555 MHZ) 和 LIb(1596MHz-1600MHz) 7j B1/L1 频段 附近 通 带 。 


L-band test data 
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Spectral channels 
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0 500 1000 1500 2000 2500 
Number of samples 
图 8 L BINE XOT 
Fig.8 L-band spectrum selection channel 
利用 上 述 所 求 阔 值 确定 遂 选 的 7 个 5MHz 无 线 电 干扰 较 少 的 无 线 电 通 带 的 洁净 率 如 下 表 4 所 示 : 
北斗 星 B1、GPSL1 频 点 北 选 通 带 为 1551MHz~1555 MHz, 1596MHz~1600MHz, 其 洁净 率 为 98.329%、 
98.301%; 
北斗 星 B2 频 点 遂 选 通 带 为 1161MHz~1165MHz、1221MHz~1225MHz, 其 洁净 率 为 98.315%、98.335%; 
北斗 星 B3 频 点 遂 选 通 带 为 1246MHz~1250MHz、1291MHz~1295MHz, 其 洁净 率 为 98.224%、97.650%; 
GPSL2 频 点 遂 选 通 带 为 1221MHz~1225MHz、1231MHz~1235MHz， 其 洁净 率 为 98.335%、98.260%， 


且 所 送 选 出 的 通 带 的 洁净 率 都 大 于 95% 以 上 ， 符 合 太 阳 观 测 需 求 。 
表 4 适 选 通 带 洁 净 率 


Tab.4 Selection belt cleanliness rate 


Cleanliness rate 


Beidou/GPS frequency point Selection pass belt 
East South West North average value 

1551MHz-1555 MHz 98.367% 98.340% 98.245% 98.363% 98.329% 

n 1596MHz-1600MHz 98.355% 98.335% 98.226% 98.286% 98.301% 
1161MHz-1165MHz 98.334% 98.194% 98.448% 98.283% 98.315% 

e 1221MHz-1225MHz 98.330% 98.351% 98.296% 98.362% 98.335% 
1246MHz-1250MHz 98.321% 98.258% 98.313% 98.004% 98.224% 

i 291MHz~1295MHz 98.21396 98.069% 97.319% 96.998% 97.650% 

L2 1231MHz~1235MHz 98.421% 98.350% 98.221% 98.046% 98.260% 


3 总 结 与 展望 


针对 即将 在 云南 天 文 台 凤凰 山本 部 建设 的 工 波段 的 精密 流量 太阳 射电 望远镜 , 本文 提出 一 种 基于 Simple 
Thresholding 算法 与 CUSUM 算法 的 改进 阐 值 算法 , 并 对 台 址 的 志波 段 无 线 电 测试 和 观测 频段 遂 选 出 介 于 北 
斗 B1、B2、B3 频 点 ,GPS L1、L2 频 点 之 间 7 个 5MHz 无 线 电 干 扰 较 少 的 无 线 电 通 带 ,分 别 为 1551MHz~1555 
MHz 、 1596MHz-1600MHz 、 1161MHz-1165MHz 、 1221MHz~1225MHz 、 1246MHz-1250MHz 、 
1291MHz~1295MHz、 1221MHz~1225MHz、1231MHz~1235MHz， 其 洁净 率 分 别 为 98.3299%、98.301%、 
98.315%、98.335% 、98.224% 、97.650%、98.335% 、98.260% ， 均 符合 太阳 观测 需求 。 为 下 一 步 接收 机 的 设 
计 和 信号 处 理 提供 了 依据 。 
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Research on observable frequency band analysis of solar radio 
current in L-band based on Simple Thresholding and CUSUM joint 


algor i thm 
YANG Hang! DONG Lang! HE Lesheng™ 


(1. Yunnan University; Kunming, —Yunnan;2. Yunnan Astronomical Observatory, Chinese Academy of Sciences, Kunming, Yunnan;) 

Abstract: L-band solar radio burst is a potential influencing factor of navigation system instability, through 
the L-band precision solar current volume monitoring can be found in real time solar burst interference navigation 
events, for this reason Yunnan Astronomical Observatory headquarters intends to establish an L-band 
multi-frequency point solar radio monitoring system. An effective assessment of the radio environment is essential 
for the stable access to observational data from the monitoring system. This paper introduces the radio monitoring 
readiness study of the monitoring platform, and conducts a 100-hour test of the L-band radio environment in the 
Fenghuangshan area of the Headquarters of the Yunnan Astronomical Observatory. In this paper, an algorithm 
based on the improved threshold of Simple Thresholding algorithm and CUSUM algorithm is proposed, and 7 radio 
passbands with less 5MHz radio interference between the Beidou B1, B2, B3 frequency points, GPS L1 and L2 
frequency points are selected, namely 1551MHz-1555 MHz, 1596MHz-1600MHz, 1161MHz-1165MHz, 
1221MHz-1225MHZz, respectively. 1246MHz-1250MHz, 1291MHz-1295MHz, 1231MHz-1235MHz, its 
cleanliness rates are 98.329906, 98.301946, 98.315926, 98.335906, 98.224906, 97.650906, 98.260906, all of which meet the 
needs of solar observation. It provides a basis for the next step of receiver design and signal processing. 


Key words: L-band; monitoring systems; radio environment; threshold; cleanliness rates 
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